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On the Anomalously Strong Natural Remanent Magnetism 
of the Lava of Mt. Utsukushi-ga-hara 


t 


By H. MATSUZAKI, K. KOBAYASHI and K. MOMOSE 
(Coll. Liberal Arts & Science, Shinshu Uniy., Matsumoto) 


Abstract 


The natural remanent magnetism of the andesitic lava of Mt. 

Utsukushi-ga-hara in the Central Part of Japan, showes that both its 
intensity and direction are too scattered to be considered as normal. 
This anomalous magnetisation is supposed to come from the mag- 
netisation owing to the momentary strong magnetic field afforded 
by thunderbolts. The writers have measured the relations between 
the demagnetisation and magnetisation of samples in a strong mag- 
netic field. The result seems to support their assumption. 

It is well known that the intensity of the natural remagnent magnetism 
(N.R.M.) of volcanic rocks is generally much larger than that of the induced mag- 
netism in the present geomagnetic field [1]. However, the lava, of which intensity 
of N.R.M. is some thirty times as large as the ordinary value, was discovered by 
_ the writers [2] on some high mountains in the Central Part of Japan as was previ- 
ously reported. 

In order to study this fact, the writers have-studied the relations between 
N.R.M. and the thermo-remanent magnetism (T.R.M.) with some pieces of rock from 
the andesitic lava of Mt. Utsukushi-ga-hara. 

At first, the samples were collected from various points close to one another 
in the same rock mass. “A-group” indicates the group of samples thus collected 
from the lower lava with an altitude of about 1400m. “B-group” those from the 
upper lava with an altitude of about 1600m. “C-group” indicates the group of 
samples collected at Ogahana, the comparatively higher point (1998m) above Mt. 
Utsukushi-ga-hara plateau. Each sample has the volume of about 20cm*. The lava 
is oliv.-bearing-hypers.-aug.-andesite which is a fairly common type of andesite in 
Japan. Though the lower lava contains much less olivine and comparatively larger 
amount of rhombic pyroxene than the upper lava. Each lava-flow shows no remark- 
able differences between them in both bulk composition and mineral assemblage as 
was noted in the previous paper [3]. 

The writers [3] measured the directions and intensities of N.R.M. and T.R.M. 
of those samples. The direction of magnetisation was measured by an astatic 
magnetomer whose magnetic moments were approximately 97 C.G.S. each. The 
sensitibity of the astatic magnetometer was calibrated by Helmholtz’ coil (72cm. in 
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Table I. Relations between N.R.M., T.R.M., Jz,¢ and Qn; of the rock samples from Mt. 
Utsukushi-ga-hara. N.R.M. of the samples (J;n), T.R.M. of the samples obtain- 
ed by cooling from Curie point, in the present geomagnetic field (J;,z) and the 
ratio between the N.R.M. and T.R.M. (Q,;) were measured. Jt. was calculated 
from J;,~- Where H=1 Oe. 


Groups No. Lacy ay Irn x a Juz x 1 bts cX rao Qni an 
511191a Ishikiriba 0.600 2.100 4.650 0.286 
511191c of 0.480 3.000 6.660 0.160 
511192 3 0.810 3.000 6.660 0.270 

“2 511193 og 0.780 4.200 9.350 0.185 
511194 Re 0.750 4.200 9.350 0.178 
511195 , ) 0.600 4.200 9.350 0.400 
519293 Hyakumagari 3.600 4.500 10.000 0.800 

B - 519295 #¢ 5.400 4.500 10.000 1.200 
519298 = 8.100 4.200 9.350 1.900 
5211206 Ogahana 97.000 3.880 8.600 25.000 
5211207 ys / 34.000 2.540 5.630 14.300 
5211208 | Fi 40.600 3.500 7.760 11.600 
5211209 - 26.200 3.020 6.700 8.650 
52112010 “5 43.500 4.500 10.020 9.600 

C 52112011 i 122.000 3.960 8.800 30.800 
52112012 3 75.000 3.940 8.740 19.000 
52112013 7; 37.800 3.560 7.900 10.600 
52112014 a 64.000 3.460 7.680 18.500 
52112015 e . 17.200 3.280 7.280 5.250 
52112016 ka 28.800 3.440 7.640 8.350 


| ¥ eS 
, BabA LEN : 
0 30° 60 dow 90° 0° Ol 6° down 90° 60° x 0° 


Fig. 1. The directions of the natural rema- 
nent magnetisation of A- and B- 
groups of samples. No. 4 indicates 
the angle of elevation, 


Fig. 2. The directions of the natural remanent 
magnetisation of C-group. No. 9, 10, 
11, 15, and 16 indicate the angles of 
elevation. 


diameter). The results are shown in Table 1 and Fig. 1. The values of both inten- 
sity and direction of N.R.M. of A- and B-group, in the same rock are scattered but 
on the whole they are the same in each group respectively. However, N.R.M. of 
the samples of C-group shows quite a random distribution of the directions, and also 
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shows values some one hundred and fifty times as large as those of A- and B-group 
concerning their intensities. For instance, the intensity of the sample No. 52112011 
(No. 11 of C-group) is 0.122e.m.u./gr. The highest value of the intensities of rock 
pieces from Mt. Ogahana that has ever been found amounts up to 0.210 e.m.u./gr. 
However, all the samples show that the intensities and directions of T.R.M. are 
nearly of the same value. It seems likely therefore that the anomalously strong 
magnetisation is more frequent in higer altitudes. 

From these facts, it was supposed that the anomalously strong values of 
N.R.M. of the samples of C-group might be intensified by such extraordinary mag- 
netic field as that afforded by thunderbolts which may be more frequent on the high 
mountains than on the lowland [1]. 

Then, the writers made an experiment to examine this assumed mechanism of 
the anomalously strong magnetisation that worked on the rocks on the top of Mt. 
Utsukusei-ga-hara. The demagnetisation of the anomalously strong N.R.M. of C- 
group by the alternating magnetic field, and the magnetisation of the same samples 
by the strong magnetic field were measured. The results show that 1) the anoma- 
lously strong N.R.M. of the samples has almost been demagnetised by the alternating 


TIRE | axygeges 
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Fig. 3. The relation between the applied 
A.C. magnetic field & the demag- 
netisation of N.R.M. (Ji-curve) & 5 
T.R.M. (Je-curve) 
O No. 21, x No. 27, @ No. 34. 


magnetic field of the intensity of 
about 400 Oe. and 2) these demagnetised 


samples have restored their intensities be 
up to the order of the anomalous 0 joo 200 300 foo 
N.R.M. by the strong magnetic field Fig. 4. Change of magnetisation with the inten- 


sities of strong magnetic field. 


of 300-400 Oe as shown in Figs. 3 and BOR ONE BONG. L 


4, And then by the same method the 
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demagnetisation of T.R.M. was also examined. The results show that T.R.M. of the 
samples is scarcely demagnetised and even the maximum value of demagnetisation 
is below 50% of the original value as shown in Fig. 3. 

From these facts the anomalously strong N.R.M. is perhaps due to isothermal 
remanent magnetism (LR.M.), so that this is due to the strong magnetisation by 
many thunderbolts which struck the same rock mass in a radom way. 

Until quite recently it has been considered that the eruption of Mt. Utsukushi- 
ga-hara took place during the Pleistocene. However, one of the writers [4] has tended 
to believe that the time of eruption could be refered to the late Tertiary or the 
earliest Pleistocene in age, based on the geomorphological evidences. 

The writers wish to express their sincere thanks to Prof. Takeshi Nagata and 
his collaborators—Mr. S. Akimoto and Mr. S. Uyeda, for their continuous and ear- 
nest encouragements during the course of this work. Acknowledgement is also 
made to the Fund of the Ministry of Education of Japan. 


References 


[1] T. Nagata, Bull. Earthq. Res. Inst., 21, pp. 2-183 (1943). 

[2] K. Kobayashi, K. Momose, and K. Kasahara, Kagaku (Science), 19, p. 284 (1949). 

[3] H. Matsuzaki, K. Kobayashi, and K. Momose, Jour. Shinshu Univ., 2, pp. 1-7. 
3, pp. 49-53. 

[4] K. Kobayashi, Geogr. Rev. Japan, 26, pp. 291-307 (1953). 


On the World-Wide Disturbance in F 2-Region 


By Tatsuzo OBAYASHI 
(Hiraiso Radio Wave Observatory, Radio Research Laboratories) 


Abstract 


The world-wide patterns of the F2-region disturbance and 
their development process with the geomagnetic storm-time were 
examined with the comprehensive world-wide data of F2-layer. 
The average disturbance in f°F2 during the main phase of the 
associated magnetic storm, so-called D(/°F2), was derived with res- 
pect to geomagnetic latitude and local-time, and it was separated 
into the Dst(f’F2) and-Ds(f°F2) component. It will be shown in 
this analysis these two components play an important role in the 
ionospheric storms, and not only their range, but also their phase 
are controlled by the magnetic activity. 
The average development course of the ionospheric storm and 
the mode of the Ds(/°F2), changing regularly with storm-time, were 
studied. The phase of the Ds(f°F2) during the active stage of 
magnetic storm is almost uniquely dependent upon local-time, but 
after the activation ceases, the pattern of the disturbance moves 
together with the rotating earth. —— 
Detailed analysis was made on the progressive aspects of the 
individual disturbances in the northern hemisphere accompanied 
with some typical severe magnetic storms. 
1. Introduction 
Many studies concerning the morphology of the average F2-layer disturbance 
have been reported by various investigators [1, 2, 3,4]. In the early stage of our 
research [5], attempts were made only to separate the time variation of disturbances 
in f°F2 and h’F2 into the two parts, the storm-time and the local-time part, using 
the data of Japanese stations. It has recently been established by D.F. Martyn [2] 
and K. Sinno [6] that the ionospheric storms at any place consist of the Ds? and 
Ds-variation (the Sp)-variation) such as defined by S. Chapman [8] in the case of 
magnetic storms. On the other hand, progressive changes of the world-wide distri- 
bution of f°F2 for individual storms also provide us with some important informa- 
tion on this subject [9, 10, 11). | 
Recently, the writer had a chance to get hourly data on ionospheric storms 
kindly sent from the various stations over the world especially at latitudes near and 


inside the northern auroral zone. 
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In this report, the world-wide distribution of the average ionospheric storm 
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and the instantaneous aspects of individual storms, will be dealt with. 


2. Ionospheric Data used in the Present Analysis 


The data examined in this study are the hourly values of /°F2 and h'F2 
during the period of magnetic storms obtained at 40 ionospheric stations scattered 
in the northern hemisphere. Ten magnetic storms are picked out here as typical 
examples, and the patterns of progressive changes of ionospheric disturbances are 
traced. The dates of these storms and some characteristics of them are given in 


subsequent Table 1. 


The abbreviations and the localities of the relevant ionospheric stations are 


Table 1. The characters of magnetic storms used in this study. 
date commence. type* maximal Kp 
Nov. 4~ 5, 1950 4th 03h ..m| ms . 6, 
LOD 5 105, =v 08 eee a ms 6, 
22~23, ,, 20 sty hoe m i 
24~25, ,, PARE 3 Sacre ms 6, 
Apr. 18~19, 1951 18, 06, 52 s £ 
Jun. 17~18, 1951 ol eee Ya s 8 
1620, 55 18, 23, 14 | ms . 6, 
Jul. Toro aie OL, eee 26 s | 8_ 
Oct. 28~29, ,, 28, . 11,754 s a. 
Dec. 27~29, ,, 27, < 21, 36 ms | 7 


for moderately severe, s for severe. 


summarized in Table 2 and shown in Fig. 1. 


Storm described by three degrees of activity; m for moderate, ms 


Table 2. List of the abtreviation and location of ionospheric stations. 


Station 


| Resolute Bay 


Point Barrow 
Tromso 
Kiruna 
Fairbanks 
Baker Lake 
Reykjavik 

| Anchorage 

| Narsarssuak 
Oslo 
Churchill 
Fort Chimo 
Fraserburgh 
Prince Rupert 
Adak 

Lindow 


' Slough 


Winnipeg 
| Fribourg 
St. John’s 


Latitude Longitude 


74.7N 
fey 
69. 
67. 
64.5 
64. 
64. 
2 
61. 
59. 
58. 
58. 
57. 
54. 
51. § 
51, 


61 


3 
9 
8 
9 
3 
1 


i) 


Station 
94.9W ) GR Giaz 
156.8,, | SC | Schwarzenburg 
19.0E OT = Ottawa 
20.555 WA) Wakkanai 
147.8W | RO | Rome 
96.0,, | AK | Akita 
218 53 WS | Washington 
149.9,, | SA | San Francisco 
45.4,, TO | Tokyo 
11L.0E WH. White Sands 
94.2W | YA | Yamagawa 
68.3,, , OK | Okinawa Is. 

ra ale TA | Formosa 
130.%5< CL | Calcutta 
176.6 ,, MU | Maui 
10.1E PU | Puerto Rico 
0.6W DK | Dakar 
97:43 TR Trinidad 
7.8E . PN | Panama 
52.7W | SG | Singapore 


Latitude 


47,.1N\ 


46.8 
45.5 
45.4 
41.9 
39.7 
38.7 
37.4 
35.7 
33.6 
31.2 
26.3 
25.0 
22.6 
20.8 
18.5 
14.6 
10.6 

9.4 

1.4 


Longitude 


(Stes 
75.8W 
141.7E 
12.5; 
140.15 
77.1W 
122.2 ,, 
139.5E 
106.5W 
130.6E 
12k. Buz9 
121553 
88.4 ,, 
156.5W 
ire ss 
t Wa: AE 
Gl2 55 
DAS ess 
103.7E 


AS.5E - 
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World Northern Hemisphere in Geographic 
Coordinates. 


On the other hand, magnetograms at 
Wakkanai Observatory, at the northern 
point of Hokkaido, are used to see detailed 
changes of the magnetic storms. The geo- 
magnetic Kp-indices are also available for 
obtaining the average magnetic activity and 
the successive storm-time changes. 
ans a 3. World-wide Aspects of the Average Iono- 


oe ee oa spheric Storm 
eat 4 In order to see the average charac- 

(oe ee teristics of ionospheric storms, the devia- 
tions of f°F2 values from their monthly 


median values during the magnetic distur- 


“9 sabe ise bances were examined. The average dis- 
a * ne coer eS turbance in /°F2, which may be called 

study. D(f°F2) variation, was derived for the 10 
typical storms at 40 stations with respect to geomagnetic latitude and local-time. 
Then it was separated into two different types of variations, the storm-time part 
Dst(f°F2) and the daily or local-time part Ds(f°F2). 

The disturbances concerned here, being denoted by Dst(f°F2), Ds(f°F2) and 
D(f°F 2), are defined as the deviations of (°F 2 during magnetically disturbed condi- 
tions while Kp indices exceed 5,, and represent the following characteristics. The 
Dst is responsible for the storm-time part of varying disturbance pattern, and 
this is essentially symmetrical about the earth’s axis. On the contrary, the Ds 
shows the effect of disturbance in modifying the average local-time variation. The 
complete pattern of disturbance corresponding to the D-variation is the combination 
of Dst and Ds components, and represents the actual systematic distortion of F2- 


Dst (FoF 2) Ds (FoF2) D (foF2) 


Fig. 2 Composition of the average ionospheric disturbance. Dst(f°F 2), Ds(f°F2) and 
D(f°F2), contours indicate the deviations of f°F2 in Mc/s. (viewed from the north 
pole) 
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layer during magnetic storms. 

These three diagrams are illustrated in Fig. 2, which may represent the 
average distributions of F2-region disturbance in the northern hemisphere (viewed 
from the north pole) during the main phase of the associated magnetic storm. In 

the above result, the following will especially be remarked. 

Dst(f°F2): The depression in Dsi is maximum near the southern edge of the 
auroral zone, and the Dsf becomes rather positive around the equator. 

Ds(f°F2): The diurnal component is predominant in Ds, its range becomes 
maximum about @=50°. The phase seems to be somewhat different with geomagne- 
tic latitude, f’F2 in the polar region being maximum about 21h~24h, and minimum 
about 9h~12h, while in middle latitudes it is maximum about 18h~21h and mini- 
mum about 7h~10h. 

D(f°F2): In the combined disturbance D=Dst+Ds, the region of depression 
in f°F2 is widely extended in comparison with that of the Ds(/’F2), but there still 
remains an area of increase in fF 2 in the dark side. 

The dependency of the D(/’F2) on season may be well represented in Fig. 3, 


D ( FoF2) 3 Summer D (foF2)* Winter 


Fig. 3 Patterns of the average ionospheric disturbance D(f’F 2) 
for northern summer and winter. 


where the left diagram shows the distribution of the D(/’F2) for the northern 
summer, and the right one for the northern winter. Remarkable discrepancies 
between these two diagrams will be seen to be that the minimum point of f°F2 is 
shifted towards earlier local-time in summer than in winter, and that the region of 
decrease in f’F2 becomes larger in summer owing to the increase in the amount 
of depression of the Dsi(/’F 2). 
4. Mode of Development of Ionospheric Storms 

The average aspect of development of the world-wide D(f°F 2) variation with 
the geomagnetic storm-time was obtained by analyzing the same data. The upper 
diagram in Fig. 4 represents the change in the magnetic activity Kp indices with 
the storm-time, averaged for the 10 storms, while the lower figures show the corre- 


sponding every 6-hourly patterns of D( f°F2) with respect to the geomagnetic lati- 
tude and local time (instead of longitude). 
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Average Storm-time Variation of Kp and D (foF2) 


6 \2 


30 36 
Storm-time 


f 12 
Fig. 4 Average storm-time changes of magnetic activity Kp and the world patterns of 
corresponding ionospheric disturbance D(f°F2), deviations of f°F2 (in Mc/s). 

As will be seen in these figures, the storm-time change in the Dst( f°F 2) which 
is independent of local-time or in other words of longitude, can be characterized 
by a slight increase in the initial phase, a remakable depression during the main 
phase and the early stage of the recovery phases, and gradual diminishing of the 
depression in the later stage of the recovery phase. The range of the diurnal 
component Ds(/°F2) is also the largest during the main phase and the early stage 
. of the recovery phase of associated magnetic storms. 

The above mentioned two facts have already been known [2, 6,7, 11]. But it 
will be noticed that the phase of the Ds(/°’F2) is delayed fairly regulary with the 
storm-time. In the early stages of magnetic storms, the phase of the ‘Ds( f°F2) is 
almost fixed to the local-time, and the depressed region of f°F2 moves around the 
earth with the sun, however, in the later stages the disturbance is not fixed to 
local-time, but is delayed gradually. In other words, the disturbance variation 
depends on local-time during the active stage of magnetic storms, but after the 
activation ceases the pattern of disturbance moves with the rotating earth, its 
magnitude gradually diminishing. It may thus be presumed that the delay of the 
phase of the Ds will be chiefly due to the transition from the active stage, where 
the Ds is subjected to local-time, to the diminishing stage of the already produced 
Ds distribution on longitude. 

5. Some Details of Individual Ionospheric Storms 

The preceding discussion deals mainly with the average characteristics of 
storms. It is also of importance to consider disturbance individually, in order to 
determine the characteristics of the irregular changes, which are eliminated when 
many storms are superposed and averaged. 

The typical storms in winter, equinox and summer seasons are selected here, 
and successive instantaneous changes of the world-wide distribution of F2-layer 
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disturbances are traced. 
1) Ionospheric Storm of December 27-29, 1951 
A moderately severe magnetic storm occurred with a sudden commencement 
at 21h 36m (U.T.) on 27th. As will be shown by the magnetic records at Wakkanai 
(; 35.7°N) in Fig. 5, after the comparatively 
long duration of positive phase in 4H, mark- 
ed decrease of H followed and the active 


state continued all day long on the 28th. 


IP XTRAS] KD 
sie a iVaPs , 
EP BIA 


A A 
24 <i Lt fghy 
gre 


21 0 3 6 9 12 15 18 21 


Dec27 ~ 28 1951 
Fic. a3 oe Afo : 
Fig. 5 Records of the magnetic storm of Dec. 27- pon on Senate piled yoisinn 
28, 1951, at Wakkanai, Japan. (time; U.T.) at 03h U.T. on Dec. 28, 1951. 


(In Figs. 6~10, 12~16 and 18~ 

20, contour lines are drawn at 

intervals of 1,0Mc/s in Af°F2) 

The ionospheric disturbance accompanying this magnetic storm was charac- 

terized by the appearance of some increased f°F2 regions localized at high latitudes 
during the early stage, and the strong development of depressed region of f°F2 in 
the later stage. The instantaneous patterns of the deviation of f°F2 from their 
monthly median values for various epochs of the magnetic storm were illustrated 
in Figs. 6 to 10. 


we 
t 


oe 2 


180° 
10 
Fig. 7 Distribution of Af°F2 during the main }}' Fig. 8 Distribution of Af°F2 during the main 
phase of magnetic storm at 09h U.T. phase of magnetic storm at 15h U.T. 


on Dec. 28, 1951. on Dec, 28, 1951 
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In the earlier phases, as will be seen in Figs. 6 and 7, positive 4/’F2 regions 
appeared near the auroral zone, grew up, and diminished intermittently or succes- 
Sively as patches at various localities in high latitudes. After these changes, 
noticeable decreased f°F2 region could be seen in the European zone. This region 
built up its activity in the sunlit side at high and middle latitudes, moving around 
the earth with the sun (Fig. 8). 

The ionospheric patterns after the declining of magnetic activity are shown 
in Figs. 9 and 10. The large depressed area of f’F2 still remained with a diameter 
of over a few thousands kilometers. Any appreciable movement of the depressed 
centre with storm-time, however, could not be found, and the disturbance was fixed 
at some particular localities. Accordingly the centre of depression was removed 
from the forenoon side to the afternoon with storm-time diminishing its magnitude 
gradually. 

. 2) Ionospheric Storm of October 28-29, 1951. 

The severe magnetic storm of October 28-29 was a typical one of compara- 

tively short duration, occurring during a rather quiet period (Fig. 11). The various 


Fig. 9 Distribution of Af’F2 during the Fig. 10 Distribution of Af°F2 during the 
last phase of magnetic storm at last phase of magnetic storm at 
Zin U.Ts on Dec. 28, 1951. 00h U.T. on Dec. 29, 1951. 


effects in the ionosphere corre- 


sponding to each stage of the 
magnetic storm were easily ob- 
tainable. 

With the advent of the 
magnetic storm, somewhat in- 


a _ 


= ane 3 2 0 3 6 9 creased region of /’F2 appeared 
Oct 28 ~ 29. 1951 


Fig. 11 Records of the magnetic storm on Oct. 28-29, at the northern American conti- 


1951, at Wakkanai, Japan. (time: U.T.) nent, as shown in Fig. 12. As 
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the magnetic storm became active, this 
positive region of 4/’F2 turned to negative 
and spread over the sunlit side near the 
auroral zone. 

During the major part of this magne- 
tic storm, the positive anomaly, which 
appeared in the west Atlantic Ocean, grew 
up rapidly and extended to the southern 
edge of the auroral zone. Concurrently, 


the negative region was also enhanced 
gradually and spread down to the lower 
latitudes. These anomalous changes, as 


= illustrated in Figs. 13 and 14, are the re- 
Fig. 12 Distribution of Af*F2 during the = markable feature of ionospheric storms at 
initial phase of magnetic storm at biotic tainted : : 
13h U.T. on Oct. 28, 1951. igh latitudes during severe magnetic ac- 


tivity. 


Wy 
i. 
5] 


ty 


Fig. 13 Distribution of Af°F2 during the Fig. 14 Distribution of Af’F2 during the 
main phase of magnetic storm at main phase of magnetic storm at 
17h U.T. on Oct. 28, 1951. 19h U.T. on Oct.-28, 1951. 


The large decrease of H with rapid fluctuations suddenly blew over at about 
20h U.T., and the gradual recovery phase followed. The positive region of 4f°F2 
diminished immediately within a few hours after the above mentioned change of 
the magnetic storm, and then a large negative region developed in the sunlit hemi- 
sphere of the temperate zone (Figs. 15 and 16). 

After the subsidence of magnetic disturbance, the distorted area gradually 
contracted, however, the weak anomalies remained for rather longer time. 

3) Ionospheric Storm of July 1-2, 1951. 

Sudden commencement followed by a major magnetic disturbance began at 
22h 28m U.T. on July 1. Moderately disturbed condition started immediately and 
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Fig. 15 Distribution of Af°’F2 during the Fig. 16 Distribution of Af°F2 during the 
; last phase of magnetic storm at last phase of magnetic storm at 
21h U.T. on Oct. 28, 1951. 00h U.T. on Oct. 29, 1951. 


active variations of the mag- 
netic field continued until about | # ‘ ~ . ~ - 
15h July 2. 

The depression of /°F2 
appeared in the night hemisphere 


soon after the outbreak of mag- 


netic storm. Then this depres- a1 0 3 6 9 2 al 21 
“a - uly 1 ~ 2, (951, 
_ sion was extended widely over Fig. 17 Records of the magnetic storm on July 1-2, 
the northern part of the world 1951, at Wakkanai, Japan (time; U.T.) 


with the progress of the storm. On the other hand, increased region of /°F2 was 
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also observed in middle latitudes at the sunset hours, so far as the earlier stage 
of the magnetic storm was concerned. Figs. 18 and 19 show this developing process 
of disturbance during the main phase of the magnetic activity. 

In the later stage, a rather diffused large negative region of 4f°F2 spread all 
over the world. While the portion of the 
greatest depression was not so clear owing 
to the blackout or masking effect by intense 
Es, the large decrease of f°F2 below f°F1 
was frequently experienced especially dur- 
ing day-time. It is also worthwhile to note 
that, in summer, such extensive depression, 
as shown in Fig. 20, remains for a consi- 
derably long time after the end of magne- 
tic storms. 

6. Conclusion 
From the results described above, it 


Fig. 20 Distribution of Af°F2 during the may be remarked that the Ds/(f°F2) and 
last phase of magnetic storm at the Ds(f°F2) variations play an important 
sal Boi ec role in the F2-region storms, and further 

that the phase of the Ds during the period of the active state of magnetic storms 

is well subjected to the local-time, but after the activation ceases the pattern of 
the disturbance moves nearly together with the rotating earth, diminishing its 
magnitude. 

Some world-wide aspects of individual ionospheric disturbances associated 
with geomagnetic storms, also confirm the above-mentioned average characters. 
Observations near the auroral zone are particularly interesting for the investigation 
of the nature of storms. Extension of this study to the investigation of the anoma- 
lous behaviour of the F2-layer near or inside the auroral zone associated with the 
intense geomagnetic variations will be desirable in the future. 
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The Effect of the Solar Tides and the Temperature Change 
on the Daily Variation in Electron Density and Height 
of the F2-Layer 


By Tatsuo SHIMAZAKI 


(Radio Research Laboratories) 


Abstract 


Chapman’s ion-production function in an isothermal ionosphere 
is first reformed for the convenience of applying it to the ionosphere 
with the temporal variation of temperature. Then using this ioniza- 
tion function, the equation of motion and the equation of continuity 
are solved, taking both the thermal- and tidal- variations into ac- 
count. As a result, the maximum electron density and its corre- 
sponding height can be obtained for each hour during a day. 
Departures of these from the norm of the static Chapman region 
indicate that the behaviour of this departure 47, is very different 
from that by Weiss [6], who deals only with tidal effect on the 
ionosphere. In our case it is inversely proportional to the tempera- 
ture variation in the course of a day except in the morning when 
the electron density is very small. It is also noted that the semv- 
diurnal variation in 4z,, appears even in the case of taking into 
account only the diurnal variation in temperature but not the tidal 
variation. As regards the behaviour of 4h,,, the result similar to 
Weiss’ can be seen, but with the reasonable difference that it is 
much raised and lowered by thermal expansion and contraction in 
the afternoon (high temperature) and in the morning (lower tempera- 
ture), respectively. Further, the effect of temperature variation and 
tidal one exerted simultaneously on the daily variations in the 
maximum electron density and its height seem to show the compli- 
cated mutual interaction between these two effects. 
1. Introduction 
Since Chapman [1] proposed a theory on the production of an ionized layer 
in the upper atmosphere, assuming an isothermal atmosphere, many attempts have 
been made to apply it to a number of ionospheric problems and with particular 
success in respect to the E and F1 regions. On the contrary, all attempts to apply 
it to the F2 region have failed to explain its complicated behaviours and, therefore, 
some improved theories have been proposed to overcome these difficulties. Among 
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these investigations, the heating effect of the sun was first introduced by Hulburt [2] 
and Appleton [3] to explain the anomalously lower electron density in summer than 
in winter and the slowly decreasing rate of electron density during the night. But 
against these considerations, Martyn and Pully [4] have asserted that the F2 region 
attains a high temperature of the order 1000°C all the year round, and Martyn [5] has 
advanced a theory of atmospheric tides, neglecting the variation in temperature. 

Later the equation of continuity for electrons subject to the influence of a vertical 
drift harmonic in time has been integrated by Weiss [6] by a numerical method. 
The recent observation of upper atmosphere by rocket measurements [7] shows 
that the systematic diurnal change in temperature in the upper atmosphere has not 
yet been firmly established, but it can naturally be expected that considerable 
daily variation in temperature may exist in the F2 region in spite of uncertainty 

‘in the lower regions. Moreover, the recent study by the writer [8] on the meridi- 
onal distribution of 4’F 2 indicates that the diurnal variation term (24 hours periodic 
term in harmonic analysis) of #’F2 is in a better distribution against the geographic 
latitude than against the geomagnetic latitude, which can be considered as the 
thermal effect in the F2 region. It is therefore worth the while that we examine 
the thermal effect on the solar ionospheric tides in the F2 region. So, in the present 
paper, the daily variation in the maximum electron density and the height at which 
this maximum occurs will be calculated by solving the differential equation by a 
numerical method, taking into account 
the diurnal variation in temperature and 
the semi-diurnal variation in vertical drift 

" velocity by horizontal tides. 

2. Derivation of the ion-production func- 
tion, taking into account the temporal 
variation in temperature 

Chapman’s ion-production function 
in an isothermal ionosphere will first be 
developed into a form fitted to the sub- 
sequent calculations in which the tempo- 
ral variation of temperature must be 
taken into account. From the preceding 
consideration, we assume as the first ap- 
proximation the atmosphere to be isother- 
mal (its temperature being T,) alliday} -eob 

ig. 


i i rtain : ' : ; P 
long in region lowers than a certa Diagram illustrating the assumed daily varia- 


height (say 200km) but to be variable in tion in temperature. It has a constant temper- 

- > ature JT. in the lower regions (4X0) during a 
temperature wl th time in the upper re- day, but has a sinusoidal variation T=T,+Ty, 
When the ’°§ -cos(wf+7) in the upper regions (h>0), being 
7 T)=T,+T,+cost at noon. + is assumed to be 
height is measured positive upwards _4,, corresponding to the minimum temperature 


from that constant height, both the at 4hr. in the morning. 


Te-Ty Te To Te +Ty 


gions, as shown in Fig. 1. 
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pressure and density remain constant at the plane of =0, because the air can 
not move across this surface by thermal expansion and contraction. The height 
distribution of air molecules at any time for the assumed atmospheric model is 
given by: 

p= pore Be , hA=<0 


= po! Few =a h>0 (2.1) 


where p is the number of molecules per cc at any height # and temperature T, 
which always equals to a constant T, where #<0 and is a function of time where 
h>0; and o is the corresponding value of p at A=0 all the time in the lower re- 
gions; and oo’ is the value at h=+0 and T=T, (initial temperature) in the upper 


regions, which equals to »)-T./T). H is “the scale height at the temperature T” - 


given by kT/mg, which is a function of only T when the constituents remain 
constant and the acceleration of gravity is assumed to be constant in the studied 
range of height. 
Then, the increase, dS, in the intensity of solar radiation incident upon the 

atmosphere (>0) corresponding to an increase in height dh is given by : 

ES To -h/ A 

dS= Apoyo ap S*Secx -e ‘/4 «dh 

Integrating this equation, we have 

S= S..+e- Adee Hosecy-e~"/ # 
So Chapman’s ion-production function can be reformed as follows: 


H h 
I(x, 4)=ABS.. “0 FyexP{ _ HW Ae0 Hosec x wa (2.2) 


which can be transformed in terms of J%' and 2’ as follows: 


I(x, h)=[e.gi-e-a™ re 
h—h™ 
Z= H x ere 


where hy is the height at which I(z,#) has the maximum value when the solar 
zenith distance is z, and J? is the corresponding value of I(x, h) at Ake. 


hy = H-In(Ap)' Hy:sec x) ; (2.5) 
m me S. 


When 2% is the corresponding value of x at {=0 (at noon, in this paper), the follow- 
ing relations can be drawn from (2.5) and (2.6): 


he Ag Sec x 

Ke. Vigwee ee Xo (2.7) 
m — Ho COS% _ 7m 
*™ H cos%z% “0 (2.8) 


2 oe 
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Then, using these relations, (2.2) is also transformed as follows: 


H, e _sec% -z 

I(x, h) =F Tye! * sects * (2.9) 
‘ a aa 
where z=2'+In eet Ty (2.10) 
When x =0, as in the case of the equator and the equinoxes, we have 
Ix, b= Pepe eee (2.11) 
' eae 

2-57 — H, (2.12) 


It is therefore concluded that the representation (2.11) similar to Chapman’s 
formula can also be adopted in our case, using the reduced height z as shown in 
(2.12). For very large values of secx (corresponding to grazing incidence of the 
beam of radiation), an appropriate correction* must be made according to Chapman’s 
theory [1]. 

3. Equation of motion and its solution 

It has been shown by Martyn [5] that the horizontal winds associated with 
atmospheric tides would produce a vertical drift motion of electrons due to the 
influence of the geomagnetic field. This vertical drift velocity v, will generally 
be a function of both height and time, but for a solar semi-diurnal drift velocity 
independent of height (umiform drift velocity) we have 


01(4, t) =V9-COS (2 wt+e) (3.1) 


where u,cose is the drift velocity at noon (/=0). 


with the diurnal variation in temperature would raise and lower the isobaric surface 
and thereby produce the vertical velocity v, of atmosphere as a whole, which 


can be represented by: 


dH _h aT 
ve(h, t= OE 7 (3.2) 


The diurnal temperature variation is now assumed to be: 


T _ T,.+T,-cos(wit+t) _ r+cos(wi+t) (3.3) 
tg T.1+T,-cost | 7y+cost ‘ 
where T= x (3.4) 


The total vertical velocity v is obtained by adding the preceding two velocities v, 


and v., that is: 


v= = oy(h, t)+ veh) (3.5) 


* Recently, corrigenda for the table of this correction appearing in Chapman’s first published 
paper [1] have been presented by himself [9] and a more accurate and detailed table proposed 


by Wilkes [10]. 
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This can be transformed by virtue of (3.1) and (3.2) 


d(h (3.6) 
Ae (2E Vo COS (2wf+€) 


‘Integrating this equation, we have 


spo = Tor; + cos <)\ coaee® dt= my HH, +cosr) [ 2-sin(wt—r+ é) 
0 


athe fig ta o7+cos(wt+r) 
—cos(2tT—€): {ror t)— 2 Da Egy (wt+r) } 
+2y-sin(2 :—6)-In{r+-cos(wt+z)} || WRT) 
Now 2 in (2.12) can be written, in terms of O(¢) in (3.7) and the abbreviation 
lyo—hy 


= , 


=O(t)+7 (3.8) 


where 7 is the reduced distance of the particular portion under consideration, 
measured from the maximum ionization level at noon, and several portions will be 
traced by following the motion, to see the daily variation in electron density and 
draw its profile against height at any hour. 

The numerical calculation of the quantity @(¢) in (3.7) reveals the non-cyclic 
variation of @(t); 0(24) does not always equal to @(0). This is due certainly to 
the fact that the raising and lowering of the ionosphere by vertical tidal drift can 
not be fully compensated with each other by thermal expansion and contraction, 
so that the initial portion of the ionospere can not return to the same level after 
the lapse of 24 hours. Under these circumstances, in order that we can hold the 
same electron density profile after 24 hours, the electron density at any height 
h will also differ from its initial value after 24 hours, corresponding to its different 
heights as mentioned above. So, practically, to find the abovementioned solution, 
we must first assume the proper electron density profile at noon so that it might 
remain constant after 24 hours under the combined influence of the temperature- 
and tidal- variations on the ionosphere during a day. But this procedure is very 
difficult to follow and it also meets with the radical difficulty that the ionosphere 
moves in one direction day by day unless there is some other mechanism than the 
assumed one, because if @(24) is larger (or smaller) than @(0), (48) is also larger 
(or smaller) than (24) on the next day and so on. 

We have thus been forced to correct the daily variation of 0 by subtracting 
a linear function, so as to make the initial value equal to the final one. As 40, the 
amount of correction required, is only about 10-20 per cent of the maximum range 
of the variation in @, the influence of this correction on the shape of the daily 
variation in ® can be neglected. Of course, this procedure must.be evidenced 
by any physical reasoning, but it may now be considered as an expedient for 


” “ - — 
a 


es a Sa A 


| 


ere 


ee 


——— 


The Effect of the Solar Tides 73 


simplifying numerical calculations. 

It is remarkable, however, that this non-cyclic variation of ® will disappear 
in the special case of cos(2r—¢)=0. But these specially limited values of the 
phase constants « and ¢ (for the diurnal temperature variation and semi-diurnal 


vertical tidal velocity, respectively) cannot reasonably be accepted as only one pos- 
sible case. 


4. Equation of continuity and its solution 


The variation in electron density N(h, ?) is controlled by the equation of conti- 


nuity : 
ON. 1-BN—div(N-2) 
dN _ ON : 
or Seg Fe: grad- ‘N=1-(8+9 Oh $i) N (4.1) 


where § is the attachment coefficient* and the ion-production function J is expressed 
as‘ in (2.11), using the reduced height z as shown in (3.8). 


Iach)=1,"- et -F,(2) (4.2) 
where angen | x'{ 90° 
=0 > kx fs90° (4.3) 


When the correction at nearly grazing incidence is necessary, sec x must be replaced 
by Chapman’s grazing incidence integral Ch (x, x) [1] [10]. By virtue of (4.2) and 
the abbreviation &,(z)=8+0v/0h, (4.1) becomes: ; 


aN 


aN Io 3 °F (2)—h,(2)-N (4.4) 


For the uniform tidal velocity, the substitution of (3.1), (3.2) and (3.5)"in k,(2) ihe 
immediately : 
hy(2)=B+ 9 SP (4.5) 
At t=0 (accordingly, T= 7, and z=7), (4.4) becomes: 
(BY) = tm Fils) hala) No) 


or Nola) = Gry: Fala) — ula) (1.6) 


where «,(7) is the special case of a,(z)= il or) (2), which is a function propor- 


tional to the time variation of N at any time ¢ and height z. Further, when 7=0, 
N,(0), the electron density at noon and at the maximum ionization level, is: 


* The removal of electrons by attachment is more probable than by recombination,.as T. Yone- 
zawa has pointed out: by his analysis of the observational data [11]. Further, when the 
attachment process, instead of the recombination one, is taken into consideration,” the differ- 
ential equation becomes linear in enon density so that we can save the labourin numerical 


calculations. 
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No 0) = Pegy Fo 0)— 2001 = (Gy L210 (4.7) 


When (4.4) is divided by (4.7): 


Oe = (0) qe gy hel) ¥ (4.8) 
where v,(2)= N.(2)/No(0) (4.9) 


From (4.6) and (4.7), the initial value of v is: 


a No(7) _ RO) F,.(7)—(%) 
= "Ne(O) ~ holm) 1—a9(0) ao 


and further, when 7=0, we have P 
vy (0)=1 (4.11) 


In carrying out numerical calculations* to solve the differential equation (4.8), 
v of the portion where 7=0 (maximum ionization level at noon) is first traced by 
following the motion with the initial condition of (4.11) and the final condition of 
having the same value after 24 hours, i.e. »2,(0)=1. This solution can be obtained 
by adjusting the quantity (0) in (4.8) to a proper value. Then, using this (0) 
and solving the equation (4.8), v of the portion at any 7 (for example 7=-+0.5, +1.0 
..) at noon is traced by following the motion under the condition that both the 
initial and final values of » should be the same as (4.10). This can be made by 
adjusting the quantity a(7) in (4.10) to a proper value. After these calculations, v 
is plotted against 7 for each hour, and from the profiles of » thus drawn we can 
obtain the maximum electron density v,, and its corresponding value 7,, for each 
hour. This 7,,(¢) is the height at noon of the portion which will have the maximum 
electron density after the lapse of time /. The actual height of maximum electron 
density at ¢ is obtained by 2,,=@(t)+7,,(¢) or more practically in units of constant 
T, by means of (3.8 ): 


ie 
Tea TAM OO+ Ff (4.12) 


Ym=Ni(Zm)/No(0) is, however, inappropriate when we intend to compare many 
cases of different assumed values of parameters with each other, because it is in 
units of N,(0), the electron density at noon and at ionization maximum level, 
possibly different in each case. So it is convenient to take the electron density in 
units of the constant /,”"/k)(0) in all cases, that is 


nie ey =? 000) (4.13) 


* The solution of (4.8). convenient for numerical calculations, can be represented by the follow- 


ing formula due originally to the method by Millington [12]. 
A= i-ale)/2_ ky(0)/2 f ‘Tor FG, © Dik Fg) | 


tat t 
1+hi(2)/2 "1 4hi(2)/2 1 Tt 1—c9(0) * Tin 1—2,(0) 


v= 


(4.7a) 
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Finally, 2,, and #,, can be obtained for each hour throughout the day. 
5. Results and discussions 

In performing the numerical analysis described in the preceding sections, 
there were used the following values of parameters: y= T,/T,=2 (corresponding to 
the ratio of the maximum to minimum in temperature of 7,+7,/7.—J,=3), r= 
%/“H,=1 (corresponding to the maximum upward drift velocity of about 25km/hr., 
assuming H)=100km), t= —4w (corresponding to the time of minimum temperature 
of 4hr. in the morning), ¢=0, 6w, 12 and 18 (corresponding to the times of 
maximum upward drift velocity of 12, 9, 6, and 3hr. respectively), /)”/H)=0.5 (cor- 
responding to the maximum ionization level at noon situated 250km high above 
the ground, assuming H,=100km), and 8=const.=0.2/hr.. 

The resulting variations in v,, and h,, are shown in Fig. 2, together with those 
in the case of y=oo and ry=0 (no thermal- and tidal- variations, or the static Chap- 
man region with the constant temperature T,) and in the case of y=2 and r=0 
(thermal variation only). It is revealed thereby that the variation in v,, is quite 
similar in all cases but that 4,, shows considerable differences of variation according 
to cases. The departures of the maximum electron density and its corresponding 
height from the norm of the static Chapman region, 4n,,=,,—7,,, and 4h,,=h,, 
—h,,,, are shown in Fig. 3 in each case of the assumed variations in temperature 
and tidal velocity. In this figure 4w,, is positive for the most part of a day and its 
behaviour differs widely from the result by Weiss [6], who deals with the case 
where only the drift motion by tides is taken into account. But these details will be 
examined in the following calculations again. The behaviour of 4h,, resembles 
that by Weiss for each phase of drift except for the reasonable difference that, in 
our case, the height is much raised (or lowered) by. thermal expansion (or contrac- 
tion) in the afternoon with higher temperature (or in the morning with lower tem- 
perature). 

In addition, the profile of electron density at noon is illustrated in Fig. 4 
together with that of the ion-production function F,(7). In each case the electron 
density is measured in units of that at the maximum ionization level which may 
differ generally with each other. It can be seen from these profiles that the height 
of maximum electron density at noon is very differently influenced in each drift 
case: it is very much raised in the cases of the drift phase (37) and (67), in 
opposition to being lowered in the cases of (97) and (127). But when the tempera- 
ture variation only is taken into account, its effect on the shape of the static 
Chapman region cannot be recognized and seems to be almost negligible : the profile 
of v in this case (y=2, r=0) is quite similar to that of the static Chapman region 
(r=™, 1a =0). 

In Fig. 5 the profile of (7) or the time derivative of the electron density at 
noon is illustrated, which shows the height of equilibrium (dN/dt=0) to be very 
high. It is to be noted that the influence of the assumed temperature variation on 
the static Chapman region at noon is to lower the time derivative of N and that 
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Fig. 3 Departures from the static Chapman region (y=oo, *=0) of the maximum 
electron density and its corresponding height, Anm=%m—Mms and Ah», 
. pare (2) 
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Fig. 4 Profiles of the electron den- 


——-f ionization Fol’. sity vs. height at noon, 
maar le o together with that of the 


(2) ionization function curve 

ae Va 12%) Fo(n) 3 vo(n) = No()/No(0) 

aa wie where 1=(M—/,”")/Ho is the 
reduced height from _ the 

maximum ionization level 

hy” at noon. 


Fig. 5 Profiles of the time derivative 
of the electron density vs. 
height at noon. 
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the tidal influences is characterized by the fact that the cases of (3+) and (6+) are 
very different from those of (9) and (12%) as may be seen in Fig’ 4, too. 

It is well known that the height of the F2 layer at middle latitudes* has two 
clear maxima at about local noon and midnight in its daily variation. This may be 
understood by assuming both of the tidal drift velocity (its phase is about (11% )) 
and temperature variation, as can be supposed by the behaviours of ,, in Fig. 2. 
The result of calculations, as illustrated in Fig. 6, indicates that the abovementioned 
supposition is correct: the behaviours of v,, and h,, in the lower two cases in Fig. 


6 seem to resemble the actual observational variations more closely than in the 


‘ ee ee 
12 15 ie 2) fe) 5, 6 9 ie 2 15 1 2 (o) 3 6 9 12 
Fig. 6 The daily variations in the maximum electron density v», and its 
corresponding height h,,; Upper—only the tidal effect is taken 
into account, its phase being (11+), Middle and Lower—both the 
tidal- and thermal- effects are taken into account. The lower 
two curves are very similar to those of the observed ones at 
middle latitudes and at the equinoxes. 


upper case (tidal variation only). Then calculating the departures of x,, and /,, 
from the norm of the static Chapman region, we can see the negative dn,, during 
the whole day in the case of (y=, v=1) taking only the tidal variation into ac- 
count, as shown in Fig. 7. This is well consistent with the result by Weiss [6]. But 


* Strictly speaking, the analysis in this paper requires y,=0, but even when 7 Na solar eee 
distance at noon) has some few value, it may be sufficient to use secy/secy) instead of sec* 
and therefore the result will need not be particularly changed. 
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Fig. 7 Departures from the static Chapman region (y=0o, v=0) of the 
maximum electron density and its corresponding height in the 
upper two cases in Fig. 6, together with that in the case of tak- 
ing only the temperature variation into account (y=2, rv=0) A 
linear superposition curve of the two cases of (y=2, r=0) and 
(y=00, r=1), in each of which only the thermal and tidal effect 
is taken into account respectively, is also illustrated. This curve 
can be seen to be very different from that of the case of (y=2, 
r=1), taking account of the two effects simultaneously. 
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when we take also the temperature variation into account, 4v,, is almost positive 
except during a short interval of high temperature in the daytime. This compara- 
tively high positive value of 4z,, is probably due to the lower temperature in the 
morening and to the fact that the rate of ion-production is very large as may be 
recognized by the relation Ie T,/T in (4.2). 

As shown in Fig. 7, 42, generally varies with time in inverse proportion to 
the temperature variation in the course of a day, except in the morning when it is 
unreasonably small in spite of the low temperature. This is probably due to the 
fact that at sunset the electron density in any case is possibly very different from 
that of the Chapman norm, but that it decreases exponentially during the night. 
Generally speaking, the difference of two exponentially decreasing functions 
(A—B)e-** becomes very small with the lapse of time, and therefore the large value 
of 4n,, is not expected in the morning. It must also be noted that the semi-diurnal 
variation in 4”,, appears even in the case of taking account only the diurnal vari- 
ation in temperature but not the tidal variation. Further, the effect of temperature 
variation on the electron density may be very favourable to explain that the electron 
density in the F2 region is depressed in the daytime, while being enlarged during 
the night, as compared with the case of the static Chapman region. 

It can also be seen in Fig. 7 that the result, by taking into account of 
both the thermal- and tidal- effects simultaneously as in the case of (7=2, v=1), is 
not equal to the linear superposition of the two cases (r=2, r=0) and (y=~, v=1), 
in each of which only the thermal- and only the tidal- effect is taken into account 
respectively. This difference may be considered as due to the mutual interaction 
between these two effects, but this must further examined in those cases of the 
temperature- and tidal- variations other than so far~considered. 

In conclusion, the writer wishes to express his cordial thanks to Mr. T. Yone- 
zawa for his suggestion and constant guidance in the course of the work and to 
Prof. K. Maeda for his special interest and encouragement. The writer is also 
indebted to Miss. T. Nameki for her assistance in carrying out the numerical 


calculations. 
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An Analysis on the Diurnal Variation of the Terrestrial Magnetism, 
especially on the day-time-variation of 
geomagnetically quiet days* 


By Masaziro OTA 
(Aso Magnetic Observatory, Kyoto University) 


Abstract 


The average feature of the geomagnetic diurnal-variation is 
conventionally described by the Fourier coefficients. Camparing these 
coeffs. of quite days and with those of all days, the characteristics 
of the S,-variation are discussed. It is recognized that at the middle 
and the low latitudes the amount of variation on quite days is large 
during the day-time and negligible during the night. If we plot a 
scheme which varies during the day-time and not during the night, 
then the average feature of the S,-variation bears a striking resem- 
blance to this scheme. The feature of this scheme is derived from 
the Fourier coeffs. computed by the conventional method, and their 
latitudinal distributions are clarified by using a peculiar latidude. 
Moreover, this scheme is very satisfactory with the potential. 
1. Introduction 
Recent studies on the geomagnetic variation_in connection with electrical 
properties of the upper atmosphere have realed that the geomagnetic variation 
observed at the earth’s surface is a synthetic effect of various properties of the 
upper atmosphere and the earth’s interior. In closely examining the results of 
world-wide observations, local effects of the observatories and unexpected anomalies 
of the S,-field are recognized. Then, if we want to discuss more delicate features 
and properties of the upper atmosphere, we must re-examine the geomagnetic data 
which have been used. For example: As Dr. Rikitake [1] has pointed out regard- 
ing the anomalous behaviour of the geomagnetic variations of short period, there 
are to be found some peculiar regions. And, Prof. Hasegawa [2] has utilized dis- 
tribution-maps of the Fourier coeffs. of the geomagnetic S,-variation as a method 
for its analysis. According to these maps, we can see several stations where 
observational data do not match that of the neibouring stations. Moreover, the 
East-component (Y) at Huancayo shows the same phase as the North-component (X) 
and a comparatively large value relative to the Magnetic Equator. This fact sug- 
gests the requirement of selecting other axes for the geomagnetic diurnal-variation. 


* This paper is a contribution of Geophysical Paper dedicated to Prof. M. Hasegawa on his 
sixtieth birth-day, January 2, 1954. 
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The writer has examined the data of the geomagnetic diurnal-variation of the 
middle and the low latitudes during the II Polar Year, using the diagrams of the 
latitudinal distribution of the Fourier coeffs. about the condition of the annual mean 
', (Summer plus Winter). In the present analysis, the geomagnetic variation-fields 
are not developped by the spherical function, but are shown by the latitudinal 
distributions of the Fourier coeffs. applyed for the local time (LT); because, we 
must calculate the terms of higher degree to converge the series. 

The present paper is aimed at the re-examination of the feature of quiet days, 
because the above mentioned facts may include the effect of a disturbing field in 
spite of the feature of quiet days. First of all, the apparent distinction between fhe 
diurnal variation of quiet days and that of disturbed days is introduced by com- 
paring with the respective Fourier coefficients of each observatory. In conclusion, 
the writer defines the feature of the geomagnetic S,-variation by a simple method. 
2. Apparent distinction between the geomagnetic diurnal-variation of quiet days 
and disturbed days 

Fig. 1 shows the observatories referred to this paper which are distributed 
within the latitudes of 60°N and 60°S. By the table of the Fourier coeffs. of diurnal 
variation of those stations (not shown in this paper), it can be seen that the first, 


Fig. 1 The Observatories distributed within 60°N-60°S, fonemiatia’ 
coordinates) 


the second and the third harmonics are predominant and that the terms higher than 
the fourth are negligible. Therefore, the feature of the diurnal variation can be 
represented by these terms, which are discussed in this paper. Figs. 2-5 show the 
latitudinal distribution of each term of the North-(X) and East-(Y) components in 
regard to quiet days and all days (and disturbed days partly). The feature of 
diurnal variation of all days seems to show the average state of disturbed days. 
The mean of diurnal variation of the International Disturbed Days is sometimes 
affected by a largely disturbed day or successive disturbance owing to a peculiar 
magnetic-storm, and then the feature of the diurnal variation is highly distorted. 
Therefore, the average feature of the diurnal variation of disturbed days is to be 
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Fig. 2 Latitudinal distributions of by; and ayy. 


dealt with the all days’ mean, but that of quiet days will be defined and represented 
by the International Quiet Days. In the diagrams of Figs. 2-7, the abscissa denotes 
the geomagnetic latitude and the ordinate denotes the force whose unit is 10~’c.g.s. 
in geographical components. The discrepancy between the geographical components 
and the geomagnetic ones of magnetic force is insignificant so far as the data of 
the middle and the low latitudes are concerned. 

The question whether the geographical, geomagnetic or other latitudes are to 
be selected in reference to the locations of the observatories is discussed in the 
following paragraph. The object of this selection is to show the same value at the 
same latitude of different longitudes. According to Hasegawa’s analysis [2] the 
geomagnetic latitude is better than the geographical one for this end. In Figs. 2-5, 
in regard to the latitudinal distributions of the coeffs., it will be noticed that the 
definite distinction between quiet days and disturbed days appears in one term in 
respect to one component, i.e., the sine-term of the first harmonic of X-comp. (x1) 
and the cosine-term of the first harmonic of Y-comp. (a@,:). About the other terms 
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Fig. 3b Latitudinal distributions of by}. 
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we cannot detect any distinctions ; as shown in Figs. 3 and 4, difference-values (all 
minus quiet) are negligible compared with their original values, or, the latitudinal 
distributions of the difference-values are rather similar to those of the quiet, though 
they should be similar to those of the disturbed. The diagram of by, and ay; of 
the disturbed days is more distinct than that of all days as shown in Fig. 2. By 
this treatment, by; and ay; can be said to he the terms of the disturbing field. Then 
we may adopt this term as a measure of the disturbance (SD). But this term takes 
a certain value even on quiet days and shows a similar diagram to that of the dis- 
turbed days as we see in Fig. 2. On the other hand, this term of quiet days 
should have a significance of the quiet state. About this contradiction, further treat- 
ment will follow in the next paragraph. 
3. The feature of the S,-variation centred at 11> in local time 

The diurnal-variation-field of the quiet days is represented by the equivalent 
current-vortex. According to Ota’s paper [3], the centre of this vortex moves around 
the earth with the sun and is situated at the middle latitude along the meridian of 
of about 115 in local time (LT). This fact means that the mode of the variation at 
the middle latitude is symmetric about 115 in LT. And, as mentioned in the pre- 
vious paragraph, the sine-term of the first harmonic of the X-comp. takes a finite 
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Fig. 4x Latitudinal distributions of @ 2 and bx. 
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Fig. 5 Latitudinal distributions of a x3 and bys. 
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value on the quiet day. If we consider a variation centred at 11h in LT. instead of 
at noon, this (asymmetric) term would almost vanish, though some value would 
remain due to error. (The situation of this centre is not necessarily on the meri- 
dian of 114 in LT.) 

As already discussed, by: can be taken as a measure of SD, the fact that dy, 
of the quiet days has a finite value is due to a phase-lag of the variation which is 
centred at noon. Figs. 6 and 7 is the latitudinal distribution of Fourier coeffs. com- 
puted through such a treatment and can be expressed as follows: 


a,,cos mt+b,,sin mt=a',,cos mt! +b’,,sin mt’ 
where : /=t415°=f+1]h (1) 
m=, 2) and 3 


As seen in Figs. 6 and 7, the sine-terms of all harmonics of the X-comp. (b'x1, D’x2 
and }’,;) and the cosine-terms of all harmonics of the Y-comp. (a’y1, @’y2. and ays) 


Fig. 6 Latitudinal distributions of a y's and b'x’s. (Fourier coeffs. of X-component 
centred at 11" LT.) 
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Fig. 7 Latitudinal distributions of a’;’s and 6’;’s (Fourier coeffs. of Y-component 
centred at 11h L.T) 


do not show any systematic feature concerning the latitudinal distribution. For the 
scattering of points on the diagrams, the writer wishes to conclude that they are 
caused by regional anomalies of the S;-field. In conclusion of this paragraph, the 
S,-variation is shown by a’,’s and B’,’s. 
4. Selection of the equator 

In the present paper, the writer does not consider the longitudinal inequalities 
of the S,-variation, and deals with the latitudinal distribution of coeffs. only. It is 
assumed that the features of the diagrams of the distributions are quite similar for 
all longitudes, and the geomagnetic coordinates are favourable for this assumption 
as the first approximation. But, precisely speaking about the latitudinal distribu- 
tion, some discrepancies are seen in the different zones; in Figs. 2-8, the black 
circles denote the stations in the European Zone, the white circles those in the 
American Zone and the crosses those in the Pacific Zone. Consequently, we wish 
to find another coordinate which is more agreeable for the different zones. Judging 
from the discrepancies between the different zone, this investigation will be resolv- 
ed by the selection of the equatorial point where the amount of variation reaches 
a decisive maximum. By Hirono’s theory [4] the range of the diurnal variation 
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takes maximum at a point where magnetic dip is equal to zero, and shows a sharp 
gradient at about 5°N & S distance from this point, as to the same meridian. 
Then, as a test-case of this investigation, the Magnetic Equator (where dip is zero) 
is adopted, i.e., the cross-point of the Magnetic Eqnhator and the geomagnetic 
meridian is taken as 0° latitude and measured by degree along the geomagnetic 
meridian. This method is also adopted by H. Maeda [5]. But this measure of 
coordinates does not satisfy the function of spherical surface, but is nearly satis- 
factory in the middle and the low latitudes. Fig. 8 are diagrams obtained by such 
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Fig. 8 Distributions of a’’s and b’y’s arranged on the modified latitude (which begins from the 
point where dip is zero) 
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a method, and we can see that the discrepancies between the different zones are 
almost ignored, compared with those of the previous paragraph (Figs. 6 & 7), and 
the feature of the latitudinal distribution at the equatorial zone is made clear in 
spite of the lack of observation at this zone. 
5. Simple expression of the S,-variation at the middle and the low latitudes 

On the analysis of the diurnal variation expressed by the harmonic coefficients, 
it is assumed that the geomagnetic force varies uniformly during a full 24 hours as 


shown in the following equations ; 


Pp 
X=@y0+)d)(@yx,,c0S mt +by,Sin ml) 
m==1 


P 

Y=a@y.+d)(ay,,cos mt+by,,sin mt) 25 
meal 

where ?¢ denotes local time (LT), a’s and b’s are functions of latitude, 
and these relations are satisfied when =3 as described in the 
paragragh 2. 


But the apparent feature recorded on magnetograms shows that the main part of 
the variations on the quiet days occurs during the day-time only. In view of this 
opinion it may be thought that the Fourier coeffs. analysed over a full 24 hours 
are not an actual existence so far as the geomagnetic S,-variation is concerned. As 
we see in Figs. 6 & 7, the feature of the latitudinal distribution of cosine-terms of 
X-comp. (a@’x1, @’x2 and a’x;) have a certain similarity, ie., these three show that 
the point of zero (where the potential is max. or min.) is situated at about 30°N and 
S on 115 meridian in LT. As it is unnatural that three elements have the same 
phases about local time and reach maximum on the same latitude, the writer wishes 
to interpret this fact as one phenomenon. The writer proposes a scheme which 
shows sine or cosine during the day-time and zero during the night, as shown by 
the thick lines in Fig. 9. Indeed, at the middle and the low latitudes, the feature 


Fig. 9a Partial cosine-curve and its harmonics Fig. 9b Partial sine-curve and its harmonics 
on a full 24 hours. on a full 24 hours. 


of X-variation on the quiet days is similar to the cosine-curve of this scheme (Fig. 
Ya), and that of Y-variation sine-curve (Fig. 9b), and they are centred at 114 in LT., 
such as shown in the equation (3); 
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X=a,'+a,'cos?’ Y= Bo’ + B,’sin r’ (3: ) 
where 7’ is defined during 55-175 and taken c’/=0 at 5b, and r’=27 at 175, 
then a)’=—a,' 8)'=0 by the assumption which X and Y are zero during the 
night (17h-5h), 
‘ and 3,’ are derived by the product of the conventional harmonic-coeffs. and the 


ay 

multiplying factor as follows; 
@ y, =0.85 a,’ b 1 = 0.44 8,’ 
@y2=—0.50,' b' 2 = —0.5 By’ (4) 
@ x3 =0.170,' b’ 3 = 0.24 By’ 


where a’s and 0’s; coeffs of the conventional harmonic-analysis 
a, and f,'; coeffs. (amplitudes) of the partial cosine- and sine-curves 


which are proposed here. 


These relations are also shown in Fig. 9. Besides, the relations are confirmed by 
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Fig. 10a Latitudinal distributions of 2's, Fig. 10b Latitudinal distributions of 84, 
> calculated by the eq. (4). calculated by the eq. (4) 
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the fact that the three harmonics of the conventional one are predominant as 
mentioned in paragraph 2. But, in this case there are three relations for a sake of 
one coeff. a,’ (or f,’), then three independent values derived from the above rela- 
tions have to coincide with each other. The latitudinal distributions of these three 
values thus obtained are shown in Fig. 10. Since by this figure the three diagrams 
which are obtained by the three relations are almost agreeable, then the question 
is solved as a whole. As a conclusion for these diagrams which are to coincide 
with each other, the writer adopts the weighted mean which is antiproportional to 
the multiplying factor of the equation (4), and it is shown in Fig. 11. 

Strictly speaking, there is some disagreement between the three diagrams. 
This disagreement is caused by terms of the full 24 hours’ waves as well as error 
in the assumption of the scheme. But as these values are of the same order of 
regional anomalies of the S,-field, the writer wishes to adopt this scheme as the 
present course of analysis. Moreover this scheme is more convenient and reliable 
than the previous one for the reason given in the next paragraph. 
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The scheme proposed here brings up two questions; one concerns the dis- 
continuities at 54 (near sunrise) and 17h (near sunset), and the other is the X-values 
during the night. About these questions the writer wishes to give the following 
interpretation; this scheme is defined during the day-time only, where the physical 
state at 5 is equal to that at 174, and the phenomenon (supposed in this paper) is 
assumed to repeat itself as a period of 12 hours. T hough the physical state during 
the night is not taken into consideration, X- and Y-components are taken to be zero 
all through the night, and to have the same values as those at 5 and 17). Then, 
to answer the questions, it must be supposed that the zonal electric current which 
is proportional to o,' (which is a function of latitude) continues throughout the day. 
6. Comparison between the coefficients derived from the X-component and from 

the Y-component 

On the determination of the coeffs. of the potential function which expresses 
the geomagnetic diurnal-variation, Prof. Champman [6] adopts the values derived 
from the Y-comp. But Hasegawa’s paper [2] expresses both values as derived from 
X and Y, the final values of which are not yet published, but which can be taken 
in the writer’s hand, as shown in Table I: 


Table I 
Main teams of the spherical harmonic representation of th” S,-field computed by X and 
Y components, after Hasegawa and Ota. (local-time-field, normalized, unit 10-*c.g.s.) 
V=RDD (A?’cos mt + BMsin mt) P70) t: local time 
mn 


cceffs. Ax Br. 

m n from X from Y from X from Y 
i 2 10.10 7.09 —2.47 —2.28 
3 3 —4,96 —4.89 1.61 2.02 
2 4 1.25 1.25 ~ 0.95 —0.78 


By this Table we can see the following results about their main terms: About 
the coeffs. of the first harmonic, those from X are larger than those from Y. About 
those of the second harmonic, both values are agreeable with each other. About 
those of the third harmonic, those from X is a little larger than those from Y. As to 
these disagreements, the writer is of the opinion that it is caused by the following 
reason; the multiplying factor for a’x, to secure o/ is 1.18 (1/0.85), and that for b’y: 
to secure f3;' is 2.23 (1/0.44) as recognized by the equation (4). Then if the present 
scheme satisfy the conditions of the potential field by which the S,-field is represent- 
ed, coeffs. of the scheme derived from X should take the same values as those 
derived from Y. About the conventional analysis the coeffs. A,'! and B.' from X are 
larger than those from Y. Though the reduction-factor does not cancel this differ- 
ence, as the writer has modified the latitude, this is almost removed. The multi- 
plying factor of the second harmonic are 2(1/0.5) for X as well as Y. This fact 
suggests that the second harmonic has a characteristic essentially the same as the 
scheme proposed here. As to the third harmonics, the writer would refrain from 
definite discussion, as their values are small compared with the first and the second. 
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The potential by the present scheme is written as follows; 


V=R(A)' + Ai‘cos 7’) 


A,! and Aj’ are functions of latitude (5) 
/ yA 
then a aoe ¥ oA cos 7! 


Y=(2/sin0)(A,/sint’) as its definition 
A,’ and at are derived from X and Y independently, using the latitudinal distri- 
bution of a,/ and f, (Fig. 11). In this calculation, values at every 5° (in latitude) 
are measured by the diagram of Fig. 11. Fig. 12 shows the latitudinal distributions 
of A,’ and REN which have two diagrams respectively. By this figure it may be 


said that the features of the two diagrams derived from X and Y almost coincide 


© Values by Y-comp. 
* Values by X-comp. 
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Fig. 12a poeraa gia re of Fig. 12b Latitudinal distributions of A’, 
< : rE “ culated by the eq. calculated by the eq. (3) and (5). 


with each other in the northern hemisphere, but not in the southern hemisphere 
where observatories are sparse and regional anomalies of the S,-field are easily 
affected in the investigation of world-wide phenomena. Therefore, the writer adopts 
the diagram of the northern hemisphere, in which mean of the values from X and 
Y (Fig. 12) is ‘used, and assumes that the field is just asymmetrical in the both 
hemispheres. Fig. 13 shows the S,-field thus obtained, and the amounts of the total 
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current by this field is about 58000 ampere. This value is very agreeable with the 
previous results. It is also investigated that this field has not a non-potential part 
by the above discussions. 


x ee Pe : al ar at P uy? 15h igh ap ain aih 
‘Fig. 13 Simplified S,-field expressd by VR(A',+A’;cos<’) © denotes the position of the 
sun. Unit: 6.4 10%c.g.s. 
7. Conclusion 

The present paper is a chapter of an analysis of the diurnal variation of the 
terrestrial magnetism and the preliminary settlement of the problem in which the 
regional anomalies, the axes of the variation and other elements related to securing 
the feature of the S,-field at the middle and the low latitudes are discussed. The 
detailed discussions about those problems will be done in the next paper which will 
be published in near future. 

The scheme of the S,-variation proposed here is not a complete expression 
in mathematics, but is able to simplify it and will be service in the investigation of 
the conductivity and the motion of the upper atmosphere. It is clear that the S,- 
field has’ no non-potential part, and it is worth examining whether the disturbing field 
have non-potential part or not. 

The terms of Z-component are perhaps important in respect to this subject, 
but their regional-anomalies are so irregular that they could not be examined. 
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Meeting of the Society of Terrestrial Magnetism and Electricity : 
The 14th General Meeting was held at the Kyoto University on Oct. 30— 
Nov. 1, 1953. 

Number of the Reports read at the Meeting: 
Cosmic Rays, 9; Atmospheric Electricity, 9; Geomagnetism, 7; Rock Magne- 
tism, 6; Radio Meteorology, 3; Ionosphere, 9. 

Tanakadate Prize was awarded for the following excellent worker: 

The 14th, Mr. K. Nagashima ; 
Theoretical Research on the Cosmic Ray Variation caused by the Electric 
Field. 


The 15th General Meeting was held at the Tokyo University on May 28— 
30, 1954. 

Number of the Reports read at the Meeting: 
Geomagnetism, 7; Ionosphere, 16; Cosmic Rays, 9; Rock Magnetism and 
Electricity, 7; Atmospheric Electricity, 14. 

Tanakadate Prize was awarded for the following excellent worker: 

The 15th, Mr. K. Hirao; 
Radio Meteological Research in the Ultra Short Wave Propagation, espe- 
cially on the Nocturnal Cooling Play. 
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